Flux reversal in a simple random walk model on a fluctuating symmetric lattice 
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A rather simple random walk model on a one-dimensional lattice is put forward. The lattice as 
a whole switches randomly between two possible states which are spatially symmetric. Both lattice 
states are identical, but translated by one site with respect to each other, and consist of infinite 
arrays of absorbing sites separated by two non-absorbing sites. Exact explicit expressions for the 
long-time velocity and the effective diffusion coefficient are obtained and discussed. In particular, it 
is shown that the direction of the steady motion can be reversed by conveniently varying the values 
of either the mean residence times in the lattice states or the transition rates to the absorbing and 
non-absorbing sites. 



PACS numbers: 05.40.Fb, 05.60.Cd, 87.16. Nn 



I. INTRODUCTION 

The study of the transport of Brownian particles along 
periodic structures has attracted considerable and ever 
growing interest in the last years due to its relevance in 
a wide variety of problems in physics, chemistry, biology, 
and nanotechnology Q, 0. As required by the second 
law of thermodynamics, Brownian motion cannot induce 
a steady directed motion in a system at equilibrium (see, 
for instance, Ref. for a detailed proof). Nevertheless, 
a Brownian particle can display a directed motion in a 
periodic potential which switches randomly or periodi- 
cally between two or more states 3, 4, 5, 6, 71. For this to 
be possible, certain spatiotemporal symmetry and super- 
symmetry conditions must be broken |8j. An important 
phenomenon often observed in these systems is the re- 
versal of the current direction upon variation of a system 
parameter 0, Hol | . Though in many cases the spatiotem- 
poral symmetry breaking is achieved by using a spatially 
asymmetric periodic potential, in Ref. |ll| it was shown 
that it can also be obtained by using a spatially sym- 
metric periodic potential which fluctuates between three 
states, one of them being of free diffusion. As shown in 
Ref. if the symmetry of glide reflection is broken, 
a directed transport can also be induced in a spatially 
symmetric potential which switches randomly between 
two states. Besides the purely theoretical interest, the 
study of this kind of models is relevant because it may 
provide valuable information about the workings of the 
so-called molecular motors p], . Molecular motors are 
able to travel along polymer filaments by utilizing chem- 
ical energy from ATP hydrolysis, and are responsible for 
the transport of different substances inside the eukaryotic 
cells. 

Recently, it has been demonstrated experimentally 
that the Brownian motion of a colloidal sphere can be 
rectified using a periodic potential that alternates be- 
tween three [lj] or two [13 states which differ only by 
a discrete translation. Each state is spatially symmetric 
and consists of a large array of optical traps, created with 
dynamic holographic optical tweezers, each of which can 
localize a colloidal sphere. The pattern of discrete opti- 



cal traps can be modeled as a periodic array of Gaussian 
wells of width a and depth Vo, uniformly separated by a 
distance L substantially larg er than a. In the three-state 
case considered in Ref. [l4(> the initial array of optical 
traps is extinguished after time T and replaced immedi- 
ately with a second array, which is displaced from the 
first by L/3. This process is repeated after time T with 
an additional displacement of L/3. A cycle is completed 
when, after time T, the array of traps is displaced once 
again by L/3 and returns to its initial state. In the two- 
state case reported in Ref. ^jj , the first step is the same 
as in |l4j |. i.e., the initial array of optical traps is extin- 
guished after time 7\ and replaced immediately with a 
second array, which is displaced from the first by L/3. 
This second array is extinguished after time T 2 and re- 
placed again by the first, thereby completing one cycle. 
One main result of the three-state case is that the direc- 
tion of the motion depends sensitively on the duration 
of the states, T, relative to the time required for the 
particle to diffuse the intertrap separation [14|. By con- 
trast, in the two-state case this direction de pen ds on the 
relative durations of the two states T 2 /Ti [T^. Specif- 
ically, the flux is directed from each optical trap in the 
longer-duration state toward the nearest optical trap in 
the short-lived state, and vanishes for T\ — T 2 . Thus, 
in both cases, by varying the value of the corresponding 
parameter the flux can be reversed. 



To gain a deeper insight into the behavior of these sys- 
tems, it would be desirable to develop simplified models, 
amenable to analytical treatment, which capture the key 
features of the original complex systems. The main fo- 
cus of this paper is to present such a simplified model 
for the experiment described in Ref. [l5|. The paper is 
organized as follows. In the following section, we intro- 
duce the model and provide definitions of the quantities 
of interest, namely, the long-time velocity and the effec- 
tive diffusion coefficient. In Sec. IIIII we obtain explicit 
expressions for these two quantities. The main results 
arising from these expressions are discussed in Sec. IIVI 
Finally, we present conclusions for the main findings of 
our work. 
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II. DESCRIPTION OF THE MODEL 

In our simplified model the motion of the colloidal 
sphere in the periodic potential is modelled by a random 
walk on a one-dimensional lattice with its sites located 
at x n = nh, with neZ, where it has been assumed that 
the distance between two consecutive sites x n+ \ —x n = h 
is independent of n. The lattice can be in two possible 
states that will be indicated by the label a = ±1. In con- 
trast with the approach described in Ref. [l5j in which 
the potential switches periodically between its two states, 
here we will assume that the lattice fluctuates randomly 
between its two possible states. The evolution of the joint 
probabilities Pn (t) that the particle is at site x n and the 
lattice in state a at time t are governed by the master 
equation 



P { n\t) 
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where 7^ is the transition rate from the lattice state a 
to —a, and Ln and Rn are the transition rates from 
site x n to site x n -i and x n +i, respectively, assuming that 
the lattice is in state a. The initial condition is irrelevant 
since we are interested in the long-time behavior of the 
system. 

Retaining the main features of the model described in 
Ref. , we will assume that both states of the lattice are 
spatially symmetric (i.e., there exists a n' £ Z such that 
L^l_„ — R^_„ for all n € Z) and periodic with period 



L ( n a) and BQ N 



R n a) for all n e Z) 



n'+n ' "n' —n 

N (i.e., L<$ N 

This is an important distinction between our model and 
other similar models previously reported in the literature 
0,0,0,0, in which at least one of the lattice states is 
spatially asymmetric. For simplicity, we will restrict our 
study to the case iV = 3, which is the smallest value of 
the period for which a non-zero steady directed motion is 
possible under the above conditions of spatial symmetry 
and periodicity [2(j. The lattice state a = +1 is obtained 
from the lattice state a = — 1 by translating it to the 
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and R 
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= R 



(-i) 



right one site, so that, L 
When the lattice is in state a = — 1, there are absorbing 
sites at the locations £3^+1, for all j G Z, and conse- 



quently, L { 3j l\ 



4jl\ = 4$ = <+2 = "■ Those 



absorbing sites play the role of the optical traps in our 
simplified model. The transition rates to an absorbing 
and a non-absorbing site will be denoted by g and k, re- 
spectively, so that, L ( ~l\ = R ( ~ 1} = L { + 1] = R^+\ = g 



and j4j+2 = Li 3j 1} = RiT } = L 3i+i = k - A sketch of 
the model just described is shown in Fig. ^ 

It is worthwhile to point out that the model just de- 
scribed is formally equivalent to another random walk 
model which mimics the motion of some molecular mo- 
tors. In this model the walker switches randomly between 
two possible states, indicated by the label a = ±1, 7W 
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FIG. 1: Sketch of the fluctuating lattice. The absorbing sites 
are depicted by crosses. 



being the transition rate from the state a to —a (see 
Fig. El . These states are characterized by the position of 
two legs, each with a head-shaped object attached at one 
of their ends. The walker "walks" on a rod which pos- 
sesses an infinite array of binding sites, depicted by grey 
squares in Fig. [21 These binding sites are separated from 
each other by a distance 3h, and translated by h/2 with 
respect to the sites of the lattice. Whenever any of the 
walker's "heads" comes in contact with a binding site, the 
walker stays trapped at that site until a fluctuation of its 
state occurs. For each state of the walker, the sites of the 
lattice with this property (absorbing sites) are depicted 
by crosses in Fig. [21 When none of the walker's "heads" 
are in contact with a binding site, the walker can hop 
to its nearest-neighbor sites. Then, the transition rate 
from a site to its neighbors is g, if after hopping any of 
the "heads" comes in contact with a binding site, or k, 
otherwise. This interpretation of the model differs from 
the one considered at the end of Ref. , in three impor- 
tant aspects. First, ours is discrete and exactly solvable. 
Second, in ours the walker fluctuates randomly between 
its two possible states, whereas in theirs it switches peri- 
odically between them. In this sense, our model may be 
more effective for describing real molecular motors, since 
in this class of molecules the directed motion arises as 
a consequence of a competition between two stochastic 
time scales (the ATP binding rate and the ATP hydroly- 
sis rate). Third, in 0] the direction of the steady motion 
depends only on the relative durations of the two states, 
whereas here, as we will see later on, it also depends on 
the sign of g — k. 

We are interested in the study of the long-time veloc- 
ity or flux v and the effective diffusion coefficient D c ^ 1 
defined as 



t^oo dt 



(2) 



and 



D eS = l lim ± [( x *(t))-(x(t)) 2 ], (3) 
2 *— >oo dt 



where the Ztli moment of the position is given by (x (t)) 



3 



State +1 




I — I X I 



X 9 I „ I .)< 



duce the quantities 
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FIG. 2: An alternative interpretation of the model sketched 
in Fig. The binding and absorbing sites are depicted by 
grey squares and crosses, respectively. 



There are certain properties of v and D c s that can 
be predicted by simple considerations. For instance, 
in the absence of fluctuations between the two lattice 
states, both the long-time velocity and the effective dif- 
fusion coefficient will be zero, as the particle ends up 
being trapped by an absorbing site. Even if the sites de- 
picted by crosses in Fig. were not absorbing and, there- 
fore, there were a non-zero transition rate from them 
(the same to the left and to the right), the long-time 
velocity would be also zero in the absence of fluctu- 
ations between the two lattice states. This is so be- 
cause, in that case, both lattice states would be in- 
variant under reflections about any of the sites with 
crosses. Another interesting property is obtained if one 
takes into account that the fluctuating periodic lattice 
sketched in Fig. ^ is invariant under the transformation 
{n, a} => {— n, — a}. Consequently, the long-time veloc- 
ity and the effective diffusion coefficient fulfill the sym- 

'-!)] = -i;[ 7 (- 1 ), 7 (+ 1 )] and 

D e g 1 7^ +1 ^l ) where we have 
written explicitly the dependence of v and Z? c ff on 7 ( +1 ) 
and 7 ( -1 ). Thus, if 7 t+ 1 ) = 7 ( _1 \ the long-time velocity 
will be zero. 



metry relations u[ 7 ' +1 '', 7 ' 
D cS [ 7 (+1) ,7 ( - 1} ] = 



and 



#°(*) = h J2 (n + l N )P ( nlw(t) ~ <z(*M a) (*). (5) 



which are periodic in n with period N, i.e., q^l N {t) 



q r \ '(t) and s^) N (t) = s%"' (t) for all n 6 Z. Substituting 
the expressions for the first and second moments of the 
position into Eqs. @ and |J3}, and replacing Pn"\t) by 
the right-hand side of Eq. |Q , it is not difficult to prove 
that the long-time velocity and the effective diffusion co- 
efficient can be expressed as 
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where = limt->oo gi a) (t) and = lim^oo s„ '(t). 
Thus, for the evaluation of v and D c r it suffices to cal- 
culate the values of the AN numbers and Sn , with 
n = 1, . . . , N and a = ±1. 

Differentiating Eqs. (QJ and © with respect to time, 
making use of Eq. as well as of the periodicity of 
and i?£ Q ' ) , and taking into account that lim^oo (t) = 
lim^oo si"" 1 (t) — 0, it is easy to see that the long-time 

limits Qn and Sn fulfill the set of AN algebraic equa- 
tions 



III. EXPLICIT EXPRESSIONS FOR THE 
LONG-TIME VELOCITY AND THE EFFECTIVE 
DIFFUSION COEFFICIENT 

In order to go further in the analysis v and -D c ff, it is 
necessary to obtain explicit expressions for them. To do 
that, we will extend the approach developed in Ref. [2l| 
to the case of a fluctuating lattice. The procedure de- 
scribed here is quite general and can be applied to any 
fluctuating periodic lattice with an arbitrary period N. 
For this reason, in what follows it will be described in 
general terms. An alternative approach based on the 
gradient expansion of the probability distribution corre- 
sponding to a coarse-grained description of the random 
walk model can be found in Ref. 17]. First, let us intro- 
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Summing up the 2N equations contained either in the 
set of equations JSJ or ©, it is easy to realize that in 
each set there are only 2N — 1 independent equations 
(excluding some pathological cases where the number 
can be even lower). Thus, to univocally determine the 
AN quantities Q„ and it is also necessary to add 

the two additional equations X)a=±i J2n=i Qn^ ~ 1 an< ^ 
J2 a =±i Zw=i ^" = 0, which can be straightforwardly 
proved from the definitions of and S^K 

To apply the above procedure to the model sketched 
in Fig. ^ first one has to solve this set of equations, with 
N = 3 and the corresponding values for the transition 
rates Ln and . Replacing the results obtained for 
Ql Q) and S ( n a) in Eqs. © and 0, one obtains after some 
lengthy simplifications that 

v = 3hA P (1 - Ap 2 ) (g - k)gkT 2 A~ l (11) 

and 

-v 2 BT~ 1 A~ 1 , (12) 

where 

A = 4g 2 {g + 2k) 2 + 8g(g + k){g + 2k)T 

+ [4(g + k) 2 - (1 + 3Ap 2 )fc 2 ] T 2 , (13) 

B = Ag 2 {g + 2k) 2 + 2Ag(g + k){g + 2k)T 

+ [28(g + k) 2 - 3(3 + Ap 2 )k 2 } T 2 + 8(g + k)T 3 , 

(14) 

r = 7 ( +1 ) + 7 (- 1 ', and Ap = pi^-pi+V, p^ 3 being the 
equilibrium population of the lattice state a, i.e., pi^ = 
ry(~ a ) jY. This quantity can alternatively be expressed as 
Ap = -T( +1 )] / [TH) where T< Q ) = 

l/<-y( Q ) i s the mean residence time of the lattice in the 
state a. 



IV. RESULTS 

An important result that follows from Eq. (|llfl is that 
the direction of the steady motion not only depends on 
the sign of Ap, as would be expected from the previ- 
ous symmetry considerations, but also on the sign of the 
difference between the transition rates to an absorbing 
and a non-absorbing site, g — k. Thus, for instance, if 
> T^ +1 \ the steady directed motion is from left 
to right for g > fc, whereas it is from right to left for 
g < k. Obviously, if < T^ +1 \ the behavior is 

just the opposite. Finally, no steady directed motion is 
possible if either = or g = k. This subtle 

dependence of the direction of the steady motion on the 



values of g and k cannot be inferred a priori from sym- 
metry considerations, and is one of the main results of 
this work. It is important to point out that this model 
predicts correctly the direction of the steady motion of 
the experiment described in Ref. ^{|. Indeed, in that 
case, the values of g and k can be estimated as the in- 
verses of mean times required for the particle to diffuse 
a length L/3 — a/2 and L/3, respectively, where L is the 
intertrap distance and a is the width of an optical trap. 
Therefore, g = 2D /(L/3 - a/2) 2 and k = 2D /{L/3) 2 , D 
being the diffusion coefficient of the colloidal sphere in 
the fluid. Since g > k, the motion is from left to right 
if r^" 1 ) > T( +1 ) and from right to left if T^ 1 ) < T( +1 \ 
i.e., from each absorbing site in the longer-duration state 
toward the nearest absorbing site in the short-lived state. 
The above expressions for g and k also show that the fi- 
nite width of the optical traps plays an important role in 
the appearance of the steady directed motion. 

Some limiting values of the expressions Hll|) and i|12|) 
can be easily interpreted. For instance, as suggested by 
the previous simple considerations, v and D e g tend to 
zero both as T — * (absence of fluctuations between the 
two lattice states) and as Ap — » ±1 (absence of fluc- 
tuations between the two lattice states in the long-time 
limit). As g — > 0, v and D e g also go to zero because, in 
this limit, the particle stays trapped between two consec- 
utive absorbing sites, since it is not allowed to jump to 
them. The expressions for v and D c g also tend to zero 
both as k — » and as g — > oo. This is due to the fact that, 
in these limits, the motion of the particle consists of a se- 
ries of forward and backward jumps from a non-absorbing 
site to its nearest absorbing site, which arise as a conse- 
quence of the fluctuations in the lattice state. Finally, 
the limiting values limr^oo v — ^oo and limr^oo -Deff = 
Ajff,oo> which can be easily calculated from Eqs. i|TT|) 
and l|12l) , also have a simple interpretation. Using the re- 
sults in Ref. [2l]], it is not difficult to check that Voo and 
-Deff,oo are the long-time velocity and the effective diffu- 
sion coefficient corresponding to a random walk with the 
average transition rates L n = piq 1 L„ ^ + pcq 1 ^! 

and R n = piq iil + piq Rn ■ In contrast with 
what happens in continuum models, the limiting value 
Voo does not vanish in general. This difference be- 
tween discrete and continuum models has been previously 
pointed out in the literature (see, for example, Ref. 
From the expressions for and -D c ff.oo, it can be eas- 
ily proved that lim^oo = 3hAp(l — Ap 2 )fc/4 and 
linig^oo £> ff,oo = 9h 2 (l — Ap 2 )k/8 which, in general, are 
not equal to zero. This reflects the fact that the lim- 
its r — > oo and g — > oo of Eqs. Ijllfl and (|12fl do not 
commute. 

In order to simplify the graphical presentation of re- 
sults, it is convenient to introduce dimensionless quan- 
tities by defining characteristic length and time scales. 
The characteristic length scale L is chosen to be the 
distance between two consecutive absorbing sites, i.e., 
L = 3h. The characteristic time scale r is chosen to 
be the mean time required for the particle to diffuse the 
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3: Dependence of the dimensionless long-time velocity 
tv (top panels) and effective diffusion coefficient L~ 2 rD e ff 
(bottom panels) on Ap (left panels) and log 10 (rg) (right pan- 
els). The dimensionless parameter values are: L~ 1 h = 1/3, 
rfc = 9/2, rg = 9 (left panels), Ap = 0.6 (right pan- 
els), tT = (solid lines), tT — 10 (dotted lines), and 
tF — 100 (long-dashed lines). The limiting values L~ 1 tv oc 
and L _2 r_D e H i00 are plotted with dash-dotted lines. 



length L in the absence of absorbing sites. Taking into 
account that in the absence of absorbing sites the effec- 
tive diffusion coefficient is -Dcff.NA = h 2 k, it follows that 
t = L 2 /(2£> cff ,NA) = 9/(2fc). Then, the number of in- 
dependent dimensionless parameters is reduced to only 
three (the values of tT, Ap, and rg), since the other two 
take the fixed values L~ x h =1/3 and rk = 9/2. 

In Fig. |3 we depict the dependence of the dimen- 
sionless long-time velocity L~ x tv (top panels) and ef- 
fective diffusion coefficient L~ 2 tD c r (bottom panels) on 
Ap (left panels) and log 10 (re?) (right panels). We have 
taken rg — 9 (left panels), Ap — 0.6 (right panels), 
and three values of tT, namely, tT = (solid lines), 10 
(dotted lines), and 100 (long-dashed lines). With dash- 
dotted lines we have plotted the limiting values L~ 1 tv oq 
and L~ 2 TD e g :OQ . Notice that both the absolute value of 
L~ 1 tv and L~ 2 tD s are increasing functions of tT, ex- 
cept at those points at which these functions vanish. This 
behavior is confirmed by a detailed study of the signs of 
the derivatives dv/dT and dD c g/dT. The reversal of the 
steady directed motion is clearly illustrated in panels (a) 
and (b). These panels also reveal that, for fixed values of 
the remaining parameters, there exist two optimal values 
of Ap (opposite to each other) and of rg which maximize 
the absolute value of the long-time velocity in each di- 



rection. In panel (b), one can also see that the location 
of the second maximum is shifted towards larger values 
of rg as tT is increased, and that this maximum com- 
pletely disappears in the limit r — > oo. Panel (c) shows 
that the dimensionless diffusion coefficient as a function 
of Ap is maximized for Ap = 0, i.e., when the mean resi- 
dence times in each lattice state are equal. Analogously, 
there is also an optimal value of rg which maximizes the 
dimensionless diffusion coefficient as a function of rg [see 
panel (d)]. The location of this last maximum exhibits 
a similar behavior with tT as the second maximum in 
panel (b). 



V. CONCLUSIONS 

In summary, we have proposed a simple random walk 
model which allows the analytical study of the flux rever- 
sal in a spatially symmetric fluctuating lattice. We have 
obtained explicit analytical expressions for the long-time 
velocity and the effective diffusion coefficient and made a 
study of them. In spite of its simplicity, by appropriately 
choosing the transition rates to an absorbing and a non- 
absorbing site, our model is able to reproduce some im- 
portant features of the experiment described in Ref . |15| . 
In particular, it provides the correct direction of the in- 
duced flux and shows that this direction is closely related 
to the finite width of the optical traps. Furthermore, our 
model predicts the behavior of the effective diffusion coef- 
ficient that would be observed in the experiment reported 
in Ref. ^||. Finally, it predicts that similar results to 
those contained in [15j would be observed if the array of 
optical traps switched randomly between its two possible 
configurations instead of periodically. 

We have also presented an alternative interpretation of 
our model which mimics the motion of some molecular 
motors. In contrast to the one considered at the end 
of Ref. pj|. we have proved that, in ours, the direction 
of the steady motion depends not only on the relative 
durations of the two states, but also on the sign of the 
difference between the transition rates to an absorbing 
and a non-absorbing site. 
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